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1 Introduction 



The agreement between LEP data |3] and standard- model predictions is overwhelming . 
Deviations from the standard model would signal new physics and are therefore intensively 
sought. For this reason, high-precision standard-model radiative corrections must be 
computed, entailing not only two-loop diagrams but also the effects of particle masses. 
In many cases, radiative corrections to leading mass-dependent effects are significant; 
in some cases [|] they can reach 50%. Thus, in recent years, the attention devoted to 
two- loop diagrams with masses has increased considerably f4]|— [fl^ . 



Whilst massless multi-loop diagrams can be dealt with by essentially algebraic meth- 
ods (IB]], the situation is more complex in the case of massive diagrams. The first successful 
algebraic approach to massive two-loop two-point functions achieved a systematic evalua- 
tion of diagrams with external massive particles on the mass-shell || [7] , using recurrence 
relations obtained from integration by parts in d = 4 — 2e spacetime dimensions. These 
recurrence relations reduced a large variety of integrals, depending on only one mass 
parameter, to combinations of T-functions and a single so-called master integral that 
contains all the problematic analytic properties of the Feynman diagrams. 

As a next step, we investigate in this paper all the master integrals, with one massive 
particle and an arbitrary external momentum, that occur in the two-loop self-energy dia- 
grams of QED and QCD, in an arbitrary dimension d. Applications of the corresponding 
d = 4 diagrams are made, for example, in QCD sum rules |14], |T5]| . There are several 
reasons for evaluating them with arbitrary d. 

First, it turns out || that <i-dimensional calculations are often much easier than their 
4-dimensional counterparts. Recurrence relations are best derived within the framework of 
dimensional regularization, where contributions from surface terms in a partial integration 
are disregarded. Moreover, these recurrence relations are systematically implementable 
as computer-algebra algorithms. Hence the algebraic methods of this paper are far less 
taxing than previous 4-dimensional analytical methods (as the author of j|] is only too 
painfully aware). 

Secondly, our methods yield not only two-loop results in QED |L6], [L7], [18|], QCD |L9|, pOp 



and QCD2 Plf , but are also applicable to the dimensional regularization of higher-loop 
calculations, where two-loop terms of order e L ~ 2 are needed to obtain an L-loop result. 

Thirdly, many of our results are analytically simpler than their specializations to 
d — 4. We find that the ci-dimensional master integrals are generally expressible in terms 
of simple hypergeometric functions; only when these are expanded in e do we encounter 
complicated polylogarithms. From this point of view, the fearsome analytic complexity 
of trilogarithms |22|], first encountered in the two- loop QED calculation of |IB| and later 



surfacing in QCD |19fl , is merely a consequence of taking a singular limit of a few otherwise 
very tractable hypergeometric functions of the type 2 -Pi and 3F2. 

Finally, the analytic continuation of hypergeometric functions giving a Taylor series in 
q 2 /m 2 , to obtain those giving an asymptotic expansion in m 2 /q 2 , is simple for arbitrary d, 
yet almost impossibly difficult if one knows only the d = 4 Taylor coefficients. Thus the 
lesser labour of determining the <i-dimensional small-g 2 expansion also yields the greater 
benefit of a pair of complementary expansions, to use for d — 4. 
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The objectives of this paper are as follows: 



(a) to obtain, whenever possible, hypergeometric representations of the <i-dimensional 
two-loop master integrals of QED and QCD; 

(b) when this appears impossible, to obtain the differential equation satisfied by the 
(^-dimensional master integral; 

(c) to generate, thereby, any desired number of terms in both the small-g 2 Taylor series 
and also the large-g 2 asymptotic expansion, for any d; 

(d) to demonstrate the utility of applying methods of accelerated convergence to each 
of the two series thereby obtained; 

(e) to exemplify these algebraic and numerical methods in the concrete case of the 
self-energy of the photon. 



These aims are complementary to those of [0 , which extends the calculation of lead- 
ing terms in the asymptotic expansion, with d — 4, to mass cases more complicated 
than those considered in ||. Here we restrict attention to the mass cases of [|J, whilst 
both significantly generalizing and also greatly simplifying the analysis, by working with 
arbitrary d. Moreover, the knowledge we derive from the complete Taylor series and 
from the complete asymptotic expansion enables us to investigate methods of accelerated 
convergence ||23|| , which we find to be of great utility. We therefore proceed as follows. 

In Section 2 we consider all the mass cases of j| , for arbitrary d. These are illustrated in 
Fig. 1. For all cases except that of Fig. 1(c), we achieve objectives (a) and (c), above. The 
irreducibility to simple hypergeometric series of integral /3 of Fig. 1(c) is to be expected; 
for d — 4 it is not reducible to polylogarithms, being rather a function whose discontinuity 
has a derivative involving an elliptic integral 0]. For J 3 , we achieve objectives (b) and (c), 
using differential equations of the type developed in El 
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Fig. 1: The mass cases of [HI: dotted lines are massless; solid lines have mass m. 







In Section 3 we study the d — 
whose Taylor series diverges for \q 2 



4 and d — > 3 limits of the most difficult integral, J3, 
> m 2 , whilst its asymptotic expansion diverges for 



\q I < 9m . A method for separating the sources of discontinuity is given in Section 3.1, 
allowing series expansions throughout the entire complex plane. This is shown to entail 
working to high numerical precision at intermediate stages of the computations. 

In Section 3.2 we find a more efficient solution to the problem of intermediate space- 
like values of q 2 , namely the application of Pade approximants to the Taylor series, and 
their extension, via the so-called e-algorithm p3[, to achieve accelerated convergence of 



the asymptotic expansion. The improvement of the convergence of the expansions, within 
their respective domains of validity, is impressive. Even more significant is the high 
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numerical accuracy achieved in the space-like region q 2 e [m 2 , 9m 2 }, where neither of the 
original series was valid. This extension of the domain of convergence is known to occur 
when one applies Pade approximants to the Taylor series of a function with a positive 
spectral density, namely to a function of the so-called Stieltjes type [^3|. We have no right 
to expect it to occur for the asymptotic expansion, since that involves both powers and 
logarithms of m 2 /q 2 , when d = 4. We regard the numerical results of Section 3.2 as highly 
significant; such dual application of methods of accelerated convergence to small-g 2 [|Il[ 
and large-g 2 series may well be the next step to take in extending the scope and 
accuracy of practical standard-model calculations. 

In Section 4 we apply the methods of previous sections to the photon self-energy. 
First we show how analytically simple are the new d- dimensional results for the self- 
energy and its spectral density. Spin is truly an inessential complication in d dimensions; 
the methods that yield the corresponding master integral also yield the full self-energy 
in terms of precisely three simple hypergeometric functions. Secondly, we show how 
impressively accurate are the results of accelerating the convergence of the Taylor series 
and the asymptotic expansion. 

In Section 5 we summarize our findings and draw conclusions from the efficiency of 
our algebraic methods and the accuracy of our numerical approximations. 



2 Master integrals 

We consider the generalization to d = 4 — 2e space-time dimensions of integrals analyzed 
in [U in the case d = 4, namely the (euclidean) integrals 

h = /J X^f+ef ^ ° )P{k2 > ° )P{kl ~ q > ° )P{k2 ~ g ' ° )P{kl ~ 0) ' (1) 
h = II X^ 2 fl+e) nK m)P{k2 > ° )P{kl ~ q > ° )P{k2 ~ q > ° )P{kl ~ h > 0) ' (2) 

h = II X^a + ef ^ 1 ' m)p ^ m)p{kl ~ q > 0)p{k2 ~ q > 0)p{kl ~ k2 > 0) ' (3) 

r c r\ tc ti Jc 

h = J J + 2 £) P(h, 0)P(fc 2 , m)P(h - q, m)P(k 2 - q, 0)P(k 1 - fc a , m), (4) 

r c f\ Ic c\ lc 

h = J J + 2 P(k h m)P(k 2 ,m)P(h- q,m)P(k 2 - q,m)P(ki- k 2 ,0), (5) 

h = II ^^^P(k 1 ,m)P(k 2 ,0)P(k 1 -q,m)P(k 2 -q,0)P(k 1 -k 2 ,m), (6) 

where P(k,m) = l/(k 2 + m 2 ) and the subscript of Jjy denotes the number of massive 
lines. The denominator structure of integrals (0-0) is illustrated in Fig. 1. (Note that the 
d — » 4 limits of these integrals, given in ||], differ by a factor of q 2 .) Integrals {I , 1 3 , 7 4 } 
are needed for the two-loop gluon self-energy; integrals {I 2 , 13} for the quark self-energy. 
As in we include the integral I\ for illustrative purposes. 
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2.1 Integral I\: an example of Mellin— Barnes techniques 



With only one massive line, the integral @) of Fig. 1(a) is most easily evaluated by using 
the Mellin-Barnes representation [ES 



P(ki, m) 



1 



k 2 + m 2 27ri 



Tfcf / dsT (- s m + s) 



m 



(7) 



where the contour separates the poles of the integrand on the right, with Res > 0, from 
those on the left, with Res < — 1. Inserting (|?|) in (|2|), one obtains an integral over s 
whose integrand is given in terms of T-functions by the triangle rule of |L3]. Closing the 
contour to the left, we find the small-g 2 expansion 



q 2 m 4£ e 3 (l-2e)I 1 = \5 2 F l 



e,2e; -g 
1-e; 



+ {13- \8) 2 F 1 



l,e; -g 2 
1-e 



+ h 2F1 



1,-e; -q 2 
1 -3e 



where 



= 

7 = 
<5 = 



r 2 (l-e) m 



2\ £ 



— = 1 -e\n(q 2 /m 2 ) + 0(e 



r(l - 2e) V g 
r(l - e)r(l + 2e) 



r(l + e) 



1 + 2((2)e 2 + 0{e 6 ) 



(9) 
(10) 
(11) 



are structures that will recur in our analysis. Closing the contour to the right, we find 
the large-g 2 expansion 



l,e; m l 



1-e 



e, 2e ; m 2 



-g 2 



-572^1 



2e,3e; m l 



i + e; -g 2 



"^7 2^1 



l,3e; m 2 



i + e; -g 2 



(12) 



which may also be obtained by inverting the hypergeometric functions of 

As d — > 4, we find that the same function appears in the small-g 2 expansion (|8]) and 
the large-g 2 expansion ([121) : 



Fi{q 2 /m 2 ) = 6C(3) + F^m 2 / q 2 ^ 



d=4 



where 



n-l 



n=l 



r=l 



n — r 



r 2 n 2 



-x 



(13) 
(14) 



for |s| < 1, in agreement with Eq. (14) of 
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2.2 Integral I4: for abelian boson self-energies 



The small-g 2 expansion of integral (|5|) of Fig. 1(d) is obtained by extending the method of 
integration by parts, used in for the massless integral ([1]). The triangle rule eliminates 
either the massless line or two of the massive lines. In the former case, one has the product 
of two one-loop integrals; in the latter case, the external momentum may be routed 
through the massless line, making the small-g 2 expansion straightforward. In terms of 
hypergeometric functions that are regular as q 2 — > 0, the <i-dimensional result is 



m 



G A 



2+4e 



e(l - 2e)h = 1 + 



4m 2 



G A — 



1,1 + e 



3 . 

2 ■ 



4m 2 



Ha 



1 'l + 2e)(l-e)' 

1,1 + e. l + 2e 



3-^2 



f + e,2-e; 
as d — ■> 4. 



4m 2 



(15) 
(16) 



which gives the small-g 2 expansion in Eq. (44) of 

This result may be transformed into one involving hypergeometric series that are 
regular as q 2 — > 00, by inverting the series (0), to obtain: 



q 2 m Ae e\l-2e)h= I 1 + 



4m 2 



-2e 



I -e, -e; 4m 
1-e; 



i\ 2 



-7 2^1 



I + s, 3e ; 4m 2 



1 + e ; -g2 



+ 2^ 2 F : 



|,2e; 4m 2 
1-e; ^g 2 



3-^2 



LI 



4m 



l-e,l-2e; -g 2 



(17) 



from which the large-g 2 expansion of Eq. (44) of 0j] is obtained in the limit — >■ 4. 



2.3 Integral I2'. for non-abelian fermion self-energies 



The key to obtaining results dT5]JT7| ) for I± was the generalization of the triangle rule of |TB|] , 
which applies when a massless particle is exchanged between particles whose masses do 
not change [pi" ]. This method is therefore applicable to integral @ of Fig. 1(b), for which 
one obtains the small-g 2 expansion 



m 2+4e e 2 (l - e)(l - 2e)I 2 = SH 2 -G 2 + , 

1 — £ m z 



with 6 given by (O) and 
G2 = 2*1 



1,1 + e; -<f 
2-e; 



#2 



2*1 



1 + e, 1 + 2e 
2-e; 



Inverting the hypergeometric series, we obtain the large-g 2 expansion 



q 2 m i£ e 3 (l-2e)I 2 =0 2 F 1 
with 



1 + e, 2e : m 2 



l-s: -a 2 



72*1 



1 + 2e, 3e ; m 
1 + e; _ 



2*1 



l,e ; m 
1-e; _ 



/3 1 + 



m 



-2c- 



1 — e m 2 



G 2 . 



(19) 



*2 + * 2 2 , (20) 



(21) 



The expansions ([0|E3) §* ve ^ ne coefficients of Eq. (22) of H, as d — > 4. 
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2.4 Integral I3: for non-abelian boson self-energies 

More difficult than either of the integrals (|3],|5]) is integral (|6]) of Fig. 1(e). In this case, 
the triangle rule fails to simplify the integral, since the exchanged particle is massive. 
Nevertheless, the small-g 2 expansion of a combination of the integral and its derivative 
w.r.t. q 2 can be reduced to 3F2 and 2F1 functions, using integration by parts. From this 
result, an expansion in terms of 4F3 and 3F2 functions is obtained by integration, which 
yields 



m 



2+4e, 



l + e)(l-2e)/ 3 



1 



'1 



2e H 



1,1,1 + e; -g 
|,2; 4m 2 



e)(l + 2e) 
2 



4-^3 



2e 



l,l + e,l + e,l + 2e; 
3 +e,2 + e,2-e; Am 2 
2 



2 1 - - 
1, 1 + e, 1 + e; -q 
§,2 + e; 4m 2 



(22) 



The presence of (m /q ) s , via (3 in the second term of ( p2|) , reveals the branchpoint at 
g 2 = 0, which prohibits the development of a simple Taylor series. Instead, one obtains 



1 



2m? 



^ T(n + ^ 



-9 



2 \ ™ 



Am 2 



n + 1 



In (g 2 /m 2 
n+1 



+ 0(e), 



(23) 



n=0 x V" 1 2; 

in agreement with Eq. (45) of 0. 

The inversion of (0), to produce a large-g 2 expansion, is likewise difficult, since each 
hypergeometric function has a repeated argument. To outflank this difficulty, one returns 
to the differential equation relating the integral and its derivative w.r.t. m? to simpler 
hypergeometric functions, whose inversion is not frustrated by repeated arguments. In- 
verting the latter and integrating the resultant series, one obtains 



q 2 m A£ e A {l - 2e)/ 3 = 4^3 
\ + £, e ; 4m 2 



1 i -£ £ 



Am 



+ 2e{l + 2e) 3 F 2 



1 + £ ; -g 2 
3 



1 — e, 1 — e, 1 — 2. ; — q 2 

1, 4m 2 
1 — £, 1 — £ ; — 



-7 3^ 



I + e, £, 3e ; Am? 



1 + 5,1 + 5; -q 2 



3-^2 



+ 



25 
1-6 



3-^2 



1,1, f + e; 4m 2 



2,2 



85 2 



4-T3 
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l,l,|,l + 2e; 4m 2 



1,1 



2,2 



1 1 2 ' 



4m 



e; -c 



2,2,2-5 



/3m 2 
g 2 



(24) 



The d — > A limit agrees with Eq. (45) of 
result of Jl3[ for integral ([[]): 



the m — > limit agrees with the massless 



g 2(d-5) F (x + gg) r 3 (-g) 

1-25 r 2 (i + 5)r(i-35) 



r(l + e)r(l-3e) 

r(i - 25) 



■ COS 7T5 



6C(3) 



+ 0(5), (25) 



which is likewise obtained from (0,0,50). 



2.5 Integral I3: for abelian fermion self-energies 

The remaining integral of Fig. 1(c) is by far the most difficult. Symptomatic of its 
intractability are the features |f] found for d — A, where the small-g 2 expansion involves 
the maximum value of Clausen's integral, which is not reducible to ^-functions, and the 
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derivative of the discontinuity in the time-like region —q 2 > 9m 2 is an elliptic integral. 
Accordingly, we generalize the method of to arbitrary d, obtaining differential equations 
whose right-hand sides do not involve intermediate states with three massive particles. 

First we treat the simplest two-point function with three massive propagators: 

d^id^ 



-P(k<z, m)P(k\ — q, m)P(k\ — m) 



7T d r 2 (l + £)" 

for which a Mellin-Barnes transform of one propagator yields 

A(n,q 2 ) + (l-2e)B(n,q 2 ) 

An 

n=0 ^ 



A{n,q 2 ) 
B(n,q 2 ) 



3 r 2 



3 r 2 



n+l,n + e,n—l + 2e; —q 
2 - e, n + \ + E] Am 2 



-1 + 2g) n 

(!+£)n ' 



n+l,n + e, n + 2— e\ —q 
3 . 

2 ' 



2-e,n + 



Am 2 



■3.) ' 

.2/" 



(26) 

(27) 

(28) 
(29) 



with (a) n = T(a + n)/T(a). At q 2 = —m 2 , we obtain from ([27H29D the on-shell result 

3 



-m A£ ~ 2 e 2 (l-e)(l-2e)J 3 



3 r 2 



l,-l+2e,i; 

\ +e,2-e; 



+ (l-2e) 3 F 2 



2 , 3 — 2e: ; 



3 

2 4" 



ie-fe 2 + 0(e 3 ), (30) 
whose e-expansion may extended, through 0(e 6 ), in terms of (Li n (l), Li n (|)|n < 5} f26fl . 



After some transformations, 



can be rewritten as 



m 4e -V(l - e)(l - 2s) J 3 = C{n) 



71=0 



(-q 2 /m 2 ) n 
(2 - e) n n\ 



C{n) 



d n 

du n 



u 



d n 

n,-l+e_~ 77i 

— 2-Ti 

du n 



l,-l+2e; u 
\ + e- 4 



+ (1 - 2e)u 



n+l-e 



du n 



l,e; it 

2 ' ^ 



(31) 



u=l 

(32) 



For the leading term in ([31]), we obtain from (p2|) the vacuum-diagram result [11 



C(0) 



7T r(2e) 



+ |(l-2e) 2 Fi 



1 

3 . 7 



2 ' 



Kl-9£ 2 5 2 ) + 0( £ 3 ), 



3 £ -l r 2 ^ 2 

where £2 @ is a multiple of Cl2(vr/3), the maximum value of Clausen's integral |22]| : 

2n- 1 4 



52 £ (3n - l) 2 (3n - 2) 2 

= 0.260 434 137 632 162 098 955 729 143 208 



t/3 / 

r- , d9 In 2 sin - 
9^ Jo V 2 



(33) 

(34) 
(35) 



and is needed to high precision, for the numerical analysis of Section 3.1. It is tedious 
to obtain higher coefficients from (|32|) . It is easier to use the second-order differential 
equation that follows from (ESLEUI), namely 



2{q 2 + m 2 )(q 2 + 9m 2 )m 2 J' 3 ' - [(d - A)q A + 10(3d - 10)g 2 m 2 + 9(5d - 16)m 4 J J' 3 



+ 3(3d -S)(d- 3)(q 2 + 3m 2 ) J 3 



(d-A) 2 ' 



(36) 



where primes denote differentiation w.r.t. m 2 . This differential equation may also be 
obtained directly, using integration by parts. It yields the recurrence relation 

C(n + 1) = l{d-2)(d-3)6 ni0 - ± {{d + A)(d-3) + Wn{d-2n-A)}C{n) 

- ±n{d + 2n-2){d-2n-2)(d-n- 2)C(n - 1) , (37) 

which enables efficient computation of any desired number of small-g 2 coefficients. 
The expansion at large q 2 is of the form 

-m 2 /q 2 ) n 

yjq un'/tj -t- pq in -fin <^ 3 yn) j — 

n=0 



q 2 m 4£ e 2 (l - 2e)J 3 = ]T (^C^n) + (3q 2 m 2 C 2 {n) + m 4 C 3 (n) 



nl 



(3* 



where Ci(0) = sj [2(1 — 3e)(2 — 3s)], C 2 (0) = 3/(1 — e), and the remaining coefficients 
are rational functions of s, determined by (^). These coefficients may be obtained in 
closed form, by transforming the hypergeometric series (|28|j29|) into large-g 2 expansions. 
This inversion generates the structures (P,|T0|). Collecting powers of —m 2 /q 2 , we obtain 
the terminating series 



Ctin) 
C 2 (n) 



;-l + 2e) n (3e) w - 2 

2(1 +£)„-! 

;-l + 2e) n (e)„ 



3-^2 



3-^2 



C 3 (n) 



(l-e)(2-e) n 
, 2(-l + 2e) B 
+ 1-s sF2 
(1 - 2 £ )n! 



1-e 



3-^2 



— ~ +e, 1 — £ — n, — n; 
-l + 2e,e; 

— I + £, — 1 + £ — n, — n : 
-l + 2£,£; 

§, 1 — £ — n, — n ; 
2 -£,£; 



5 -£, 1 -£ - n, -n; 4 
2 - e, 3 - 2e ; 



(39) 



(40) 



(2-£)„ 3 2 



|, — 1 + £ — n, — n : 
2 - £,£; 



(41) 



for the coefficients of ([$1 
equation (|36|) . 

From the expansions of J3, we generate those of 
d^d^a 



which we have checked by substitution in the differential 



K« = 



7l d T 2 (l +£ 

using integration by parts, which yields 



P(fci, 0)P(fc 2 , m)P(A; 2 - q, m)P(k 1 - k 2 , m), 



(42) 



(q 2 + m 2 )K 3 - \{q 2 + 9m 2 )4 + (3d - 8) J 3 



Am 



2d- 10 



(d-3)(d-A) \ d-A 



m 2 (q 2 — m 2 ) ^ (q 2 + m 2 ) 2 G 2 



d-2 



(43) 



with G 2 given by (|i~9|) at small g 2 and by ( f2lf ) at large q 2 . Finally we obtain I 3 from 
integration by parts, which yields 



(d - A)((q 2 + 3m 2 )I 3 - \(q 2 - 3m 2 )/ 2 ) - m 2 (q 2 + m 2 )(2I 3 + l' 2 ) + K 3 - f 3 



Am 



2d-12 



2m 4 ^ m 2 (q 2 + m 2 )G 2 q 2 (q 2 + m 2 )G 2S 



d-3 Ud-A) 2 (d— 2)(d— 4) 



(d-2)= 



(44) 



with I2 given by ( |T8"| , |2T)D . From this first-order equation the expansion coefficients may be 
obtained, at both large and small q 2 , for any d. 

The small-g 2 expansion is of the form 



m 



4e+2 r .2 



5 ~ (-q 2 /2m 2 ) n ~ (q 2 /18m 2 ) n 



e 2 {l-e){l-2e)h = 5 £ V ?' J Pi(n) + (3(0) £ J P 2 (n 



2 n ^o (d/2)„n! iV 7 w ^ (d/2)„n! 

-StS^-UL-ts^' (45) 

where 5 is given by ([Tl|), C(0) is the constant of (p3|), the final product is set to unity for 
n — 0, and P;(n) are polynomials in cZ. Table 1 gives the first 5 terms in the Taylor series. 
(We have computed the first 20.) 



Table 1: Expansion coefficients of 7 3 at small q 2 . 

Pi(0) = l, Pi(l) = d 2 -9c/ + 22, P x (2) = c/ 4 - 21c/ 3 + 168c/ 2 - 604c/ + 832, 
P x (3) = c/ 6 - 36c/ 5 + 541c/ 4 - 4338c/ 3 + 19588c/ 2 - 47304c/ + 48000, P x (4) = d s - 54c/ 7 
+ 1271c/ 6 - 17022c/ 5 + 141884c/ 4 - 754056c/ 3 + 2497696c/ 2 - 4722432c! + 3916800, 

p 2 (0) = 1, P 2 (l) = d 2 + 7c/- 66, P 2 (2) = c/ 4 + 15c/ 3 -28c/ 2 - 2484c/ + 10944, 
P 2 (3) = c/ 6 + 24c/ 5 + 133c/ 4 - 5574c/ 3 - 62540c/ 2 + 1058280c/ - 3248640, P 2 (4) = c/ 8 + 34c/ 7 
+ 439c/ 6 - 6854c/ 5 - 278612c/ 4 + 309976c/ 3 + 69223392c/ 2 - 615536640c/ + 1505018880, 

P 3 (0) = |, P 3 (l) = d 3 -7d 2 + lQd-6, P 3 (2) = c/ 5 + 16c/ 4 -463c/ 3 +2960c/ 2 -7308c/+7488, 
P 3 (3) = c/ 8 +25c/ 7 +193c/ 6 -24269c/ 5 +313846c/ 4 -1598060c/ 3 +3568632c/ 2 -2875392c/-839808, 
P 3 (4) = c/ 10 + 35c/ 9 + 509c/ 8 - 4471c/ 7 - 1276222c/ 6 + 28901060c/ 5 - 252836984c/ 4 
+ 1041981408c/ 3 - 1924952832c/ 2 + 880201728c/ + 1209323520. 

The large-g 2 expansion is of the form 

00 / 2 / ^2\n 

qVe 2 {\ - 2e)h = £ ((3 2 M (n) - 7 Mi(n) + /3M 2 (n) + M 3 (n)) ^ (46) 

n=0 n - 

where the coefficients are given by the recurrence relations 

(n + e)M (n) = n(n - 1 + 3e)M (n - 1) + (-1 + 4e) n , (47) 
(n + e)M 1 (n) = 2n(n - 1 + 2e)M l {n - 1) - n(n - \){n - 2 + 3e)M 1 (n - 2) 
+ |Ci(n + 1) + (n - 2 + 3e)Ci(n) 

+ |(2n - 1 + 3e)(-l + 2e) n (3e) n _i/(l + e) n , (48) 
nM 2 (n) = 2n(n-l + e)M 2 (n-l)-n(n-l)(n-2 + 2e)M 2 (n-2) 
+ §(n + 1 - e)C 2 (n) + n(n - 2 + 2 £ )C 2 (n - 1) 
+ (2n - 1 + e)(-l + 2e) n (e) n _ 1 /(l - e) n 



-2(-2 + 2 £ )„, 3 P 2 



l,e,-n; 1 
1 — e, 3 — 2e — n ; 



(49) 



(n-e)M 3 (n) = 2n(n - l)M 3 (n - 1) - n(n - l)(n - 2 + e)M 3 (n - 2) + 2M 4 (n) 

- §n(n + 1 - 2e)C 3 (n - 1) - n(n - l)(n - 2 + £)C 3 (n - 2), (50) 
(ra-2e)M 4 (n) = n(n - 3 + 2e)(M 4 (n - 1) - 5 n>2 ) 

- 2e(l - 2e)(n! - 5„,o - |M( £ W(2 - e) n . 1} (51) 
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with M 0i i(0) = 1/s, M 2 , 3 ,4(0) = 0, and Ci, 2)3 (n) given by (gggggg). The extra set of 
coefficients in ( |5lD serves to expand the square of series (|21|), needed for the 7 2 term on the 
left-hand side of (44]) and for the G2 term on the right. Relations (f47| --|51~1) allow very rapid 
computation of the asymptotic expansion (f4"6|); the hard work has been done in obtaining 

) of the irreducibly difficult integral fl26|) , 
). Table 2 gives the first 5 terms in the 



the closed forms (39-41), for the expansion 
together with the differential equations 
asymptotic expansion. (We have computed the first 20.) 



Table 2: Expansion coefficients of 7 3 at large g 2 . 



M (0) 
M (3) 
M (4) 

Mi(0) 
M x (3) 
Mi (4) 

M 2 (0) 
M 2 (3) 
M 2 (4) 

M 3 (0) 
M 3 (3) 



M (l) 



4(d-5) 
(d-6) 



M (2) 



-4(2d 3 -33d 2 +176d-308) 
(d-6)(d-8) 



(d-4) 

8(2rf 3 -37d 2 +210d-384)(d-5)(d-7) 

(d-6)(d-8)(d-10) ' 
-8(4d 5 -144d 4 +1987d 3 - 13158d 2 +42064d-52416) (d-5) 
(d-8)(d-10)(d-12) ' 

-2 Mi(l) — 12 ( d ~ 4 )( d ~ 5 ) />/f. — -6(5d 3 -86d 2 +472d-800) (3d- 14) (d-5) 



M x (2) 



(d-4)' '"A-V — (d-6) 2 ' ' U1 ^ ~~ (d-6) 2 (d-8) 2 ' 

12(7d 4 -195d 3 +1952d 2 -8172d+11712)(3d-14)(3d-16)(d-5)(d-7) 

(d-6) 2 (d-8) 2 (d-10) 2 ' 
-18(41d 6 -2046d 5 +41532d 4 -436680d 3 +2489536d 2 -7211904d+8128512)(3d-14)(3d-16)(d-5)(d-7) 

(d-6)(d-8) 2 (d-10) 2 (d-12) 2 

2(d-3)(d-8) A/f / \ _ -2(d 3 -15d 2 +22d+152)(d-3) 

4(2d 5 -45d 4 +197d 3 +108d 2 +3500d-23040)(d-3) 

3(d+2)(d-2)d ' 
-4(2rf 7 -51d 6 +90d 5 +3063d 4 -6180d 3 +42420d 2 - 1300496d+4660992) (d-3) 
3(d+4)(d+2)(d-2)d ' 

-8(d 2 -19d+44) 



0, M 2 (l) 



-24(39rf 2 -482d+1072) 
(d+2)(d-2) 2 d ' 



M 3 (2) 
M 3 (4) = 



(d-2) 2 d 

192(d 5 +21d 4 -824d 3 +7104d 2 -11672d-8352) 
(d+4)(d+2)(d-2) 2 d 2 



Taking the limit d — > 4, we obtain, from Tables 1 and 2, the expansions given in 



for 



q 2 h 



d=4 



00 4wdw 



w 2 — m 2 



\ w w 2 — m 2 , w 2 — m 2 \ w 2 + q 2 
In — In — In — 



+ 



4dw 



2m (w 2 — 4m 2 ) 1 / 2 



3 In 



m 



w 



w- 



3w 2 -3m 2 + g 2 W+ + W. 



m 



W+W- 



In 



w, - w. 



In 



w 

2m ' 



(52) 



with W± = ((w ± m) 2 + q 2 ) 1 / 2 . (Note that the sign of q 2 in Q] is opposite to that chosen 
here.) In the case d — 3, we find that J 3 is given by the imaginary part of a complex 
dilogarithm, throughout the space-like region q 2 > 0, and is hence reducible to three 



instances of Clausen's integral [22 



4 T 



2ImLi 2 (re w ) _ 2culnr + Cl 2 (2w) + Cl 2 (20) - C1 2 (2cj + 20) 



d=3 



rtan(0/2) 



rtan(0/2) 



(53) 



with 



m 2 + g 2 
4m 2 



cos0 



2 2 

m — g 
m 2 + g 2 



cos a; 



3m 2 + g 2 



(9m 2 + g 2 ) 1 / 2 (m 2 + g 2 ) 1 / 2 



(54) 
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3 The problem of intermediate q 2 



In the cases of integrals Ji, J 4 , J 2 and J3, the small-g 2 expansions (|,|l^,|T^,^) and the 
large-g 2 expansions ([T^Jr^j20|j2^) together cover the entire complex g 2 -plane, cut along 



the negative real axis, with the possible exception of the circle \q 2 \ = m 2 , in the cases of 
I\ and J 2 , or the circle \q 2 \ = 4m 2 , in the cases of J 4 and I3. On these circles, the series 
may converge only conditionally, or may diverge, depending on the value of d. Otherwise, 
one has a method for computing the integrals by truncation of the appropriate series, for 
any d and q 2 . 

This happy circumstance is unfortunately far from generic. In the general mass case, 
the integral 

h = J J ^dpin _|_ 2 £ ) - P ^ 1 ' m i)- P ( /c 2, m 2 )P(k 1 - q, m 3 )P(k 2 - q, m 4 )P(ki - k 2 , m 5 ), (55) 

has four contributions, with distinct branchpoints, corresponding to the four ways of 
cutting the master diagram jS. The small-g 2 expansion converges for 



\q I < Min ((mi + m 3 ) , (m 2 + m 4 ) , (mi + m 4 + m 5 ) , (m 2 + m 3 + m 5 ) J , (56) 
whilst the large-g 2 expansion converges for 

\q 2 \ > Max ((mi + m 3 ) 2 , (m 2 + m 4 ) 2 , (m 4 + m 4 + m 5 ) 2 , (m 2 + m 3 + m 5 ) 2 ) . (57) 

In the case of a branchpoint at the origin, as occurs for Ii and ^3, a small-g 2 expansion 
may be found, valid up to the next branchpoint, as in (PU2"2"|). In the general mass case, 
however, each method fails at intermediate values of \q 2 \. In particular, J3 cannot be 
computed directly from the results of Tables 1 and 2 when q 2 is in the problematic 
annulus 9m 2 > \q 2 \ > m 2 . 

Tables 3 and 4 show the accuracies to which truncations of the small and large q 2 series 
approximate the exact results (|52]j53|) , and their analytic continuations to the time-like 



region — q 2 G [0,m 2 ]. They show the relative error, A^(n) = |/3 runc //3 — 1|, when J^ rnnc is 
obtained from the first n terms in the series, with d = 4 or d = 3. 

The domain of convergence of the small-g 2 Taylor series of Table 3 is \q 2 \ < m 2 , for 
d = 4, and \q 2 \ < m 2 , for d = 3. There is a branchpoint at q 2 = — m 2 , where a cut 



begins, due to intermediate states with only one massive particle. The limit of (|52|) exists 
as q 2 -> -m 2 , giving § m 2 J 3 -> 6((2) In 2 - §C(3), for d = 4. The limit of (||) does not 
exist; instead one finds that m 2 (m 2 + q 2 )I 3 — > 4^(2), for d = 3. Note that the convergence 
of the d = 4 Taylor series is slow for \q 2 \ = m 2 . 

One expects the asymptotic truncation errors to become substantial for q 2 w 9m 2 , 
since the intermediate states with three massive particles generate a branchpoint at q 2 = 
—9m 2 . The last row of Table 4 shows that one may trespass only a small distance inside the 
annulus horribilis before the divergence of the asymptotic expansion becomes a problem. 
The task of computing the integral reliably, by series expansions at intermediate q 2 , is 
clearly a delicate one, which we accomplish in Section 3.1, by separating the sources of 
discontinuity. In Section 3.2, we show how to avoid this separation, for space-like q 2 , by 
using methods of accelerated convergence [^3[ . 
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Table 3: Small-g 2 truncation errors. 



T 



A 4 (10) 

2.3 xlO" 2 
2.6 xl0~ 3 
1.5 xlO~ 4 
2.8xl0~ 6 
2.9xl0~ 9 

3.4 xl(T 9 
4.0 xlO- 6 
2.7xl0~ 4 
6.4xl0~ 3 
1.7X10" 1 



A 4 (15) 
1.2xl0- 2 
4.4xl0~ 4 
6.1xl0~ 6 
1.5xl0~ 8 
4.8xl0~ 13 
5.8xl0~ 13 
2.2 xlO- 8 
1.2xl0~ 5 
1.2 xl0~ 3 
1.2X10' 1 



q • jm 



A 4 (5) 



6.7xl0" 2 
2.3xl0- 2 
5.6xl0~ 3 
7.7xl0- 4 
2.5xl0~ 5 
2.9xl0~ 5 
l.OxlO- 3 
8.8 xl0~ 3 
4.6xl0- 2 
2. 7x10^ 



A 3 (5) A 3 (10) A 3 (15) 



1.0 
0.8 
0.6 
0.4 
0.2 
-0.2 
-0.4 
-0.6 
-0.8 
-1.0 



5-lxlO- 1 1.7X10- 1 5.7xl0~ 2 

1.2X10- 1 9.3xl0- 3 7.3xl0~ 4 

1.5xl0~ 2 1.6xl0~ 4 1.6xl0~ 6 

4.5 xlO" 4 1.5 xl0~ 7 4.8 x 10~ n 

4.2xlO" 4 1.4xl0" 7 4.4xlO~ n 

1.3xl0- 2 1.3xlO~ 4 1.4xlO- 6 

9.1xl0~ 2 7.2xl0~ 3 5.7xl0~ 4 

3.6X10- 1 1.2X10" 1 3.9xl0~ 2 



Table 4: Large-g 2 truncation errors. 



q 2 /m 2 


A 4 (5) A 4 (10) A 4 (15) 


A 3 (5) A 3 (10) A 3 (15) 


15 
13 
11 
9 
7 


1.8xl0~ 5 2.5 xlO" 7 6.4xl0~ 9 
4.8xl0~ 5 1.3xl0^ 6 6.7xl0~ 8 
1.5xl0~ 4 8.6xl0" 6 l.OxlO- 6 
5.3xl0~ 4 8.1xl0~ 5 2.5xl0" 5 
2.6 xl0~ 3 1.3 xlO^ 3 1.4xl0~ 3 


1.8xl0" 3 3.3xl0~ 5 l.lxl0~ 6 
3.8xl0" 3 1.4xl0~ 4 l.OxlO" 5 
9.1xl0~ 3 8.0xl0- 4 1.3xl0" 4 
2.6xl0~ 2 6.2xl0~ 3 2.7xl0~ 3 
9.9xl0" 2 8.1xl0~ 2 1.2X10- 1 



3.1 Separating the discontinuities 



Denoting ||] the two contributions to ( |52|) as I a and lb, we find separate asymptotic 
expansions of the form 



I a , b = 3C(3) =F - 2C(2)L + 7C(3)) + £ ( a -^L 2 + b -^L + Ca > b{n) 



2 1 

-m 



i v n n" n 3 J \ q 2 



(58) 

with L = \n(q 2 /m 2 ). The coefficients of the simpler term, I a , with no contribution from 
intermediate states with three massive particles, may be obtained in closed form: 



n-1 n/ n 



a a (n) = -§, b a (n) = -3, c a (n) = -4 - 2n 2 C(2) + £ - (2 - -J , (59) 

from which we obtain the coefficients of h by subtraction, given the d — > 4 limit of the 
asymptotic expansion ([46]). The small-g 2 expansions are of the form 

^[ n \m 2 J ~ym \ r ) 



from which we obtain the coefficients of lb from the d — > 4 limit of the Taylor series (45). 

The strategy now is clear: in the annulus one uses the large-g 2 expansion of I a and 
the small-g 2 expansion of (The small-g 2 expansion of I a and the large-g 2 expansion of 
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Ib are never needed; outside the annulus one deals with the full integral.) A computa- 
tional difficulty becomes apparent, however, when one evaluates the small-g 2 coefficients 
of the troublesome term, If,: the coefficient db{n) is numerically undetermined unless the 
constant fl3~j|) is known to an accuracy substantially better than 1 part in 9 n . For example, 



4(20) = -2C(2) 



57479895945135907790704519109851 _|_ 934497714974357881 



14345652399118288857618205996800 



33354363399333138 



6.074 X 10 



-22 



(61) 



is 22 orders of magnitude smaller than its ^-dependent term, 
relative contribution to the Taylor series is 2 x 10~ 4 . 



Yet, at q 2 = 9m 2 , its 



It may now be seen that the method of [11], which gives the Taylor series of Feynman 



integrals in terms of S*2, requires one to work with intermediate numerical precision much 
greater than the desired final accuracy: 18 orders of magnitude greater in the above 
example. Fortunately one may exploit this novel feature and use the near vanishing of 
higher coefficients in the series to obtain a suitably accurate approximation to 5*2, which 
is how we computed the value (^), using the analytical expression for (4(30). By high 
precision computation of the appropriate series, we then obtained the acceptably small 
truncation errors of Table 5. 



q 2 /m 2 



Table 5: Intermediate-g 2 truncation errors. 



A 4 (5) A 4 (10) A 4 (15) 



l.OxlO" 3 8.7x10- 



1.5xl0" 5 
2.0 xl0~ 6 



A 3 (5) . 

7.6 xlO" 3 1.3 xlO" 3 3.4 xlO" 4 

3.6xl0~ 3 3.0xl0~ 4 4.1xl0~ 5 

2.0xl0- 3 5.7xl0~ 5 3.6xl0~ 6 

2.1 xl0~ 3 8.4xl0~ 6 2.0xl0~ 7 

5.9xl0~ 3 6.9xl0~ 6 1.2xl0~ 8 

3.0xHT 2 1.4xl0~ 4 5.9xl0~ 7 

3.2X10" 1 l.lxlQ- 2 3.6xl0~ 4 



As(10) A 3 (15) 



7 
6 
5 
4 
3 
2 



4.9 xl0~ 4 2.2x10 

2.2 xl0~ 4 4.5xl0~ 6 1.9xl0~ 7 

9.9xl0~ 5 7.0 xl0~ 7 1.2xl0~ 8 

8.8xl0~ 5 1.0xl0~ 7 4.2 xl0~ 10 

3.6 xlO" 4 7.8xl0~ 7 2.2xl0~ 9 

4.0 xl0~ 3 6.9 xl0~ 5 1.5 xl0~ 6 



For d = 3, one must separate (^) as follows: 

m A T r , r r Amr + 4C1 2 (7T - g) - C(2)tan(g/2) 
mI * d=z -G a + G b , G a - rtan (0/2) (62) 

and use the asymptotic expansion of G a and the Taylor series of G^ in the annulus. For 
the latter, high precision is again required, to deal with cancellations of the type 

A/n\ , 9291381626641020913518612265142 _ 1656090499861696 I 3 ^ i nyr? y i n-22 / fi o\ 
^ v / 3909280494259534670432360629425 166966608033225 2 ~ ' 1 

which occurs in the coefficient analogous to 



3.2 Accelerated convergence 

Another strategy for dealing with the problematic region q 2 e [m 2 , 9m 2 ] is the use of 
techniques that accelerate the convergence of the large and small q 2 expansions. A conve- 
nient technique is the e-algorithm of P^fl . In general, given a sequence of approximations, 
{S n \n = 0, 1, 2, . . .}, one constructs a table of approximants using: 

T(m, n) = T(m -2,n + l) + 1/ {T(m - 1, n + 1) - T(m - 1, n)} , (64) 
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with T(0, n) = and T(— l,n) = 0. If the sequence {S n } is obtained by successive 
truncation of a Taylor series, the approximant T(2k,j) is identical to the [k + j/k] Pade 
approximant f23|, derived from the first 2k + j + l terms in the Taylor series. Moreover, in 
the case of Stieltjes series, with a constant-sign spectral density, the n — > oo limit of diag- 
onal Pade approximants converges throughout the space-like region, giving the analytic 
continuation of the function whose derivatives at the origin are the Taylor coefficients. 
We have verified that 



)2L53) are of Stieltjes type. 



Table 6: Accelerated convergence of the Taylor series of ^3, with d = 4. 



q 2 /m 2 


£4(2) 




£4(4) 




£4(6) 




£4(8) 




3.0 


6.0x10" 


-3 


5.5x10" 


-5 


5.8x10" 


-7 


5.9x10" 


-9 


2.0 


1.9x10" 


-3 


7.5x10" 


-6 


3.4x10" 


-8 


1.5x10" 


-10 


1.0 


1.9x10" 


-4 


1.3x10" 


-7 


1.0x10" 


-10 


8.2x10" 


-14 


0.8 


8.3x10" 


-5 


3.0x10" 


-8 


1.3x10" 


-11 


5.2x10" 


-15 


0.6 


2.7x10" 


-5 


4.1x10" 


-9 


7.0x10" 


-13 


1.2x10" 


-16 


0.4 


5.0x10" 


-6 


2.0x10" 


-10 


9.3x10" 


-15 


4.3x10" 


-19 


0.2 


2.3x10" 


-7 


8.2x10" 


-13 


3.2x10" 


-18 


1.3x10" 


-23 


-0.2 


5.7x10" 


-7 


4.7x10" 


-12 


4.2x10" 


-17 


3.9x10" 


-22 


-0.4 


3.2x10" 


-5 


7.3x10" 


-9 


1.8x10" 


-12 


4.4x10" 


-16 


-0.6 


5.1x10" 


-1 


1.1x10" 


-6 


2.8x10" 


-9 


6.5x10" 


-12 


-0.8 


5.7x10" 


-3 


1.1x10" 


-4 


2.2x10" 


-6 


3.8x10" 


-8 


-1.0 


1.2x10" 


-1 


4.6x10" 


-2 


1.8x10" 


-2 


4.3x10" 


-3 



Table 6 shows the errors E±{n) = I/3//3 — 1|, with J| obtained, via the e-algorithm, 
from T(2n, 0), i.e. from the [n/n] Pade approximant to the first 2n + 1 terms of the Taylor 
series of J3, with d — 4. Comparison of Tables 3 and 6 reveals the accelerated convergence 
for \q 2 \ < m? . Table 6 also shows the extension of the domain of convergence to q 2 > m 2 . 

For larger values of q 2 , we apply the e-algorithm to the asymptotic expansion of 1%. 
Now the nth term in the series is (—m 2 /q 2 ) n times a quadratic polynomial in \n(q 2 /m 2 ), 
and hence there is no Pade approximant corresponding to the first 2n + 1 terms, nor any 
theorem, known to us, to guarantee convergence of the e-algorithm for \q 2 \ < 9m 2 , outside 
the domain of convergence of the asymptotic expansion. Nevertheless, remarkably good 
results are obtained in the space-like region q 2 > 3m 2 , as can be seen from Table 7. 

Comparison of Tables 4 and 7 reveals accelerated convergence, in an enlarged domain, 
as was observed at small q 2 . Together, Tables 6 and 7 show that results accurate to 
a few parts in 10 9 can be obtained throughout the space-like region, q 2 > 0. Thus, if 
one needs only space-like values, the delicate separation method of Section 3.1 may be 
avoided; it suffices to use the e-algorithm, to achieve accelerated convergence of the large- 
q 2 expansion, for q 2 > 3m 2 , or the small-g 2 expansion, for q 2 < 3m 2 . Only if one requires 
the real and imaginary parts on the problematic region of the cut, with — q 2 = s + iO and 
s G [m 2 ,9m 2 ], does one need to disentangle the separate sources of discontinuity. 
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Table 7: Accelerated convergence of the asymptotic expansion of I3, with d = 4. 



q 2 /m 2 


£4(2) 




£ 4 (4) 




£4(6) 




£4(8) 




12.0 


8.9x10" 


-7 


2.2x10" 


-9 


5.7x10" 


-14 


9.0x10" 


-17 


11.0 


2.0x10" 


-6 


7.6x10" 


-9 


8.6x10" 


-14 


5.7x10" 


-16 


10.0 


4.2x10" 


-6 


4.5x10" 


-8 


8.0x10" 


-13 


6.0x10" 


-15 


9.0 


9.0x10" 


-6 


1.4x10" 


-7 


4.1x10" 


-12 


3.6x10" 


-14 


8.0 


2.0x10" 


-5 


9.7x10" 


-8 


2.0x10" 


-11 


4.6x10" 


-14 


7.0 


4.6x10" 


-5 


1.4x10" 


-7 


1.1x10" 


-10 


1.7x10" 


-13 


6.0 


1.2x10" 


-4 


2.7x10" 


-7 


6.7x10" 


-10 


4.2x10" 


-15 


5.0 


3.3x10" 


-4 


6.1x10" 


-7 


5.3x10" 


-9 


1.4x10" 


-12 


4.0 


1.1x10" 


-3 


1.4x10" 


-6 


6.1x10" 


-8 


4.6x10" 


-11 


3.0 


5.3x10" 


-3 


1.8x10" 


-6 


1.7x10" 


-6 


1.7x10" 


-9 



4 Photon self-energy 

Evaluating two-loop diagrams for the photon self-energy, one encounters the integrals 
d d k 1 d d k 2 



P n (k 1 ,m)P n (k 2 ,m)P ri (k 1 - q,m)P' J4 (k 2 - g,m)P J5 (&i - fc 2 ,0) 

(65) 



7r d T 2 (l + e 

= v Q {ju 32, 33, k, h) i"' 2 )'' ' ' ?i 1 ' 



with the propagators of Fig. 1(d) raised to integer powers ji to j'5. Integration by parts 
yields the recurrence relations 



T(j 5 ,j2,ji)V Q 
T{js,ji,k)V Q 



{j 2 2 + (5" 
{jil + (5" 



1- )+j 4 4 + (5--3-)}^ 

2- )+j 3 3 + (5--4-)}Vq 



(66) 
(67) 



where T(j k ,j h j m ) = d - 2j k - ji - j m and l ± V Q (j 1 , ...) = V Q (ji ± 1, ...), etc. Using 
to remove one propagator, carrying the momentum ki, and then (|67|) to remove another, 
carrying the momentum k 2 , every Vq integral can be reduced to combinations of a rational 
function of q 2 /m 2 and d (deriving from tadpole diagrams, through which q does not flow) 
and rational multiples of integrals of the type 



2, 1,0,0) 
1^(2,2,1,1,0) 

Vg(0,a,/9,0, 7 ) 



Ga 



4-cT 

4 Lt 4 , 



v{a + (3 + 1 - d)r(f - 7)rpg + 7 - f )r(a + 7- 
r( a )r(/5)r(f )r( a + /? + 27 - d)r 2 (3 - §) 



4-^3 



7,a + /5 + 7-rf, /? + 7-f,« + 7 

^ ~ _l g+P-d ™ 1 «+/3-rf+i . 
2 ' ' ' 2 ' ' ~ r 2 ' 



d . 
2 ' 



4m 2 



(68) 
(69) 



(70) 



Moreover, the 4 F 3 hypergeometric functions of ([70]), with integer values of a, (3 and 7, 
are always reducible to two contiguous 3 F 2 functions. 



It follows that every two-loop calculation of the perturbative JL9], ^(J contribution, or 



of a gluon-condensate contribution, to the correlator of a current connecting quarks of 
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equal mass, amounts to no more than the evaluation of 5 rational functions, corresponding 
to the coefficients of unity, G4, G 2 , and two suitably chosen 3 F 2 functions. In the case 
of the master integral J4, with unit exponents, the functions G4 and H± were chosen 
in ([16]) in such a way that only G\ and H4 appeared in the result (|T5|). In general, a basis 
{l,Gi,Gl, 7 4 , J4} gives coefficients that are regular as d — > 4, enabling one to obtain as 
many terms as desired in the small-g 2 expansion, or the large-g 2 expansion, merely by 
taking the d — > 4 limits of the 5 rational coefficients, in this basis, and using the d = 4 
expansions of G4 and J4. In this sense, J 4 is truly the master integral, for this mass case. 



4.1 Vacuum polarization function 



In the on-shell renormalization scheme of conventional QED [16|, |17], |18| , the renormalized 
photon propagator has a denominator (1 + 11), where the vacuum polarization function, 
II, vanishes at q 2 = 0. Generalizing to arbitrary d, we write 



U(z) 



T(d) 



a 



(71) 

where z = —q 2 /Am 2 , traces are normalized by 1^(7^7^) = T(d)g^ u , and we form the 
dimensionless coupling a = T(3 — d/2)m d_4 e 2 /(47r) d / 2_1 from the on-shell charge and 
mass. The renormalization-scheme-invariant coupling is given by «/(l + II). For d = 4, 
a = e 2 /47r is the fine structure constant of QED, measured at q 2 = 0. 

To give a compact form for ci-dimensional two-loop vacuum polarization, we choose 
the basis {1, <px, <pf, 2 , 03 }, where 



Z Z 2-^1 



1,4 



2 ' 7 
2 ' 



(d - 3)(d - A)(d - Q) 
(d + 2)d(d-5)(d-7) 

(d-Q) 
(d + 2)d(d- l)(d-7) 



3-^2 
3-^2 



l,4-|,5-d: 

_i_ 9— d . 

2 ' 2 ' 

1,4 — |,6 — d; 

o 1 a 9—d . 
2 ' 2 ' 



(72) 
(73) 
(74) 



We implemented (|6q-|70|) in both Reduce, running on a VAX4100, and Form, 
running on a AT486. In each few seconds of CPUtime yielded 



n 2 (*) 
+ 



{z(d-2)-[l + z(d-2)]&(z)} 



(75) 



3(d — 1, 

-32z f rf 4 - 12rf 3 + 48d 2 - 53d - 24 z(d - 2)(d 2 -7d+ 16) [1 - 4>i{z)f 



9(d-4) 



d(d-l)(d-5) 



[d-6- (d-3)&(;z)] 5 



(rf-2)(d-5)(rf-6) 0iW _(^ 0?w 



4(d-l) 
2 ^2 



2(d-3)(l-z) 2(d-l)(d-3) v y 2(d-l) 

9 ^ i 8 [(d-2)(rf-5)-^ 2 -5d + 8)] 1 

2(2;) + 1 — ^ M z ) ?■ 



d-3 



d-3 



(76) 



The d — > 4 limit of n 2 (,2) is best taken after transforming to the basis {1, G4, G 2 , J4, J4}, 
in which the coefficients are regular, though rather lengthy (on account of factors of (3c?— 8) 
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and (3c? — 10) in their denominators). At d = 4, with T(4) = 4, we obtain 



20 ± _ 4(1 -,)(! + 2z) 
9 3z 3z v ; ' 



n 2 (2) 



d=4 



10 26 4(l-z)(3 + 2z) 

3 32 2; 
2(1 + 22) " 



G(z) + 



2(1 - z){\ - I62) 
3~2 



G 2 (^) 



3z 



l + 2z(l-z) 



dz 



Hz) 



G{z) 

Hz) 



G 



4ld=4 



1,1; 
3 . 

2 ' 



2ylny 



(77) 

(78) 
(79) 



q 2 h 



d=A 



6 [C(3) + 4 LiaC-2/) + 2 Li 3 (y)] - 8 [2 Li 2 (-y) + Li 2 (y)] In?/ 



2[21n(l + 2/ ) + ln(l- 2/ )]ln 2 y. 



l-l/2-l 



I/2 + I 



(80) 
(81) 



with polylogarithms []52| in ( p0[ ) that were obtained in [T^j and later derived dispersively 
in [f|]. In the next section we derive the discontinuities of ( ff8| , |S0"D from our hypergeometric 
functions 

For d = 3, 112(2) has an infrared divergence, attributable to on-shell singularities in 
the fermion propagator, which prohibit on-shell mass renormalization. A finite result may, 
however, be obtained in terms of the bare mass. 



For d = 2, with T(2) = 2, we obtain 

16 2(1 - 22) 



4 2 2 . . 
~-~ + ~G(z) 
3 2 2 



n, 2 



cZ=2 



3 4 :(l - :) 



G(z) 



2(1 - 2) 



G 2 (z) . 



(82) 
(83) 



4.2 Spectral density 

For the vacuum polarization function 11(2), we write the spectral representation 

p(t)dt 



U(z) = 2 



1 t — 2 



with spectral density 



p(t) = Imn^+W) = I&O { ( « ) Pl(t) + f « ) ft(t) + 0(a 3 ) 



Tit 

From the general integral representations 

Tic) 



Ait 



1F1 



3-^2 



a, 6 ; 

c; 

a, 6, c ; 
g?, e ; 



(t - iy- b ~H a - c dt 

(t - z) a 



r(6)r(c-6) ji 

r (d) r (e) 
r (6) r ( c ) r (d + e - b - c ) 

'A (t - 2) 



e — c, e — 6 ; 
d + e — 6 — c : 



l-t 



54) 



*5) 



(86) 



57) 
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we obtain the spectral representations of the hypergeometric functions in ( [72] , |73| ,|73]) , which 
all have a = 1. To obtain the contribution to the spectral density of the <j)\(z) terms 
in (|76|) , we use the transformation 



iF 1 



a, b ; 
c; 



r (c) r (c - a - b) 
T(c-a)T(c-b) 



a, b ; 
1 + a + b — c : 



1 - z 



+ z 1 - c (l-z) 



_ a _ b T(c)T(a + b-c) 



1 — a, 1 — b ; 
1 — a — b + c ; 



r (a) r (6) 

which enables one to split <pi(t + iO) into its real and imaginary parts, for t G [1, oo]. In 
this particular case, with a = 1, the second hypergeometric function in (^) is trivial. 

Our results for the (i-dimensional spectral density are as follows: 

l + t{d- 2 



Pi(t) 

P2(t) 



2t 2 



9) 



2* 2 pi W - r 2 (t) - ^— r 3 (t) 



(rf-3)(rf-4)t 2 [ " riv ^ d- 1 

2 + t(d- 2){d - 4)(d - 5) - t 2 (d - 2)(d 2 - 7rf + 16) 



n(t) 

r 2 (*) 
r 3 (*) 
r 4 (*) 







2(d 


r G) 


t~3(t - 1)T- 






fi)r(3-|) ' 




nt 


r(^)H(t 


- !) 2 


r(3 


d)r( M - 5 )r 




r (C 


r(^)H(t 


. 3d-9 

- !) 2 


r(5 


rf)r( M ~ 7 )r 


(3-1) 



2-^1 
2-^1 



r (g) r (g) sin (I ) 



ri(t) + 2^1 



d— 1 1 . 
2 ' 2 ' 1 _ f 

2 ' 



d-3 1 . 
2 ' 2 ' 1 _ f 

2 ' 

1,3 

7-d . 

2 ' 



2 ' 1 -t 



ri(t)r 4 (t)j, (90) 
(91) 

(92) 

(93) 

(94) 



where 7*2,3 (t) derive from the imaginary parts of (^1,0), whilst 7*4,1 (i) derive from the real 
and imaginary parts of (0), with the real part entering via lm(0 2 ) = 2Re(0i) lm(0i). 

We now investigate the limits as d — > 2 and d — > 4. The former case was considered 



in the QCD 2 sum rules of the latter was the subject of the classic QED calculation 



of Kallen and Sabry [16|. 
As d — > 2, we obtain 



2-7)3' 



r 4 (t) 



2tv 2 At 2 v 3 



1 — v 1 — 



with 

The spectral density is obtained, from ( |39"|J)0"| ) , as 

t 3 P (t) 



(95) 
(96) 









d=2 


i 





l + v 2 (1-v 2 ) 2 N 
-2^ + ^^ ln ^ 
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(£:) 2 + °(« 3 )- ( 97 ) 



Making the transformations z — > t + iO, y — > —x + iO, Iny — > \nx + in, in ([79]), we verify 
that (^) corresponds to the discontinuity of the d = 2 result of (^,^3|). Transforming 
to the bare charge and mass, we verify the QCD2 result of [^TJ. Since p(t) has a non- 
integrable singularity at t = 1, the dispersion relation should be written for (z — 1)11(2;). 
The corresponding inclusion of an additional factor of (t — 1) in the spectral integral will 
not destroy its convergence at large t. 



To obtain the d = 4 spectral density of fL6[ , we need to expand the hypergeometric 
functions of (|92]-|94]) to first order in e = 2 — d/2. The required expansions are 



d-l 1 . 

2 ' 2 ' 1 + 
3d-5 . 1 ~ 1 

2 ' 



15^/x 

7 T^xY- 



(l_ x 2 )( I +x+ 12 ) 



- 4e In 



1 + x 



+ (x + x 2 + x 3 ) [lnx + ef(x)] - e(^-x + %x 2 + fx 3 + x 4 ) lnx + 0(e 2 ) 



,(98) 



2^1 



d-3 1 . 

2 ' 2 ' 1 + 
3d-7 . 1 — 1 

2 ' 



— 3a/x 
2(1 -x\ 



X 



1 - 3e - 4e In 



1 + x 



+ (1 + x 2 ) [In 2 + ef{x)\ - 2^(3 + x + 5x 2 ) lnx + 0(s 2 
-2x 



1 ^ — - ■ 

2 ' 

/(x) = 61nxln(l • 



1 



;i + 2e)lna;-ep(x) + 0(V 



(99) 



(100) 



x) 



±ln 2 x + 6Li 2 (x) -4C(2) + 4Li 2 (-x) + 41n21nx, (101) 



g(x) = 2 In 2 ln(l — x) — ~ In 2 x + 2 Li 2 (x) 
from which we obtain 



2C(2), 



(102) 



tp(t)\ 



2v(3-v 2 ) 



a 



^ ) + { 2,(5 - 3^ 



+ 



16w(3 — v 2 



1+v 2 d 

h x— 

2v ax 



[l-v 2 )(15 + v 2 
3 



\nx 



+ 0(a 3 ), (103) 



in terms of the dilogarithmic discontinuity 

L 2 (X) : 



J 4 (_g2 / 4m 2 = t + iQ ) _ J 4 (_ g 2 / 4m 2 = t _ iQ ) 

m t 

2vri 

4Li 2 (x) + 2Li 2 (-x) + [21n(l - x) + ln(l + x)\ In 2. 



d=4 



(104) 



which was found by explicit integration in |1| and has here been obtained from the O(e) 
terms of the hypergeometric functions (|98|- |100| ). Differentiating (|104|) in our compact 



result (|103|) , we easily reproduce the lengthier formulae of [16], [17| [18| . 

Here we have obtained the d = 4 result of ( |103| , |104| ) by developing the ^-expansion 
of the hypergeometric functions in our algebraically computed d- dimensional result (|90|). 
This is the generic method we propose for unsolved problems, such as the two-loop gluon- 
condensate contributions to equal-mass current correlators ||1 5|| . 

Note that our use of the regular basis {1, G4, G\, J 4 , 1'^}, to obtain (|78|) in the previous 
section, reduces the d = 4 spectral problem of this section to the calculation of the discon- 
tinuity of the master integral (|I~5|), for d = 4. Once this is obtained from hypergeometric 
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functions, as in ( |104[ ), all other spectral calculations of this type are reduced to algebra. 
Once one has used Cauchy's theorem to construct the function with this discontinuity, 
as was done in || in the case of (|30D, all calculations of similar two-point functions are 
reduced to c/-dimensional algebra, amounting to no more than using recurrence relations 
to find 5 rational coefficients that are regular as d — > 4. 



4.3 Accelerated convergence of the Taylor series 

For small q 2 /m 2 = —4z, we find the coefficients of the Taylor series 



U 2 (z) = J2 c(n) [-Az] n , (105) 



n=l 



from the hypergeometric series of ([72| , |73| , [74]) that enter our d- dimensional result (|76[) . 



Running Reduce 3.4.1 on a Vax 4100, it took 20 seconds of CPUtime to obtain the 
first 10 Taylor coefficients, the numerator of c(10) being a polynomial in d of degree 27. 
We list only 

4 (3c/ 5 - 48c/ 4 + 227c/ 3 - 58c/ 2 - 1728c/ + 1872) 
C{ } ~ 9(c/ + 2)c/(c/-3)(c/-5)(c/-7) ' ( } 

lie/ 7 - 227c/ 6 + 1289c/ 5 + 1619c/ 4 - 31840c/ 3 + 30036c/ 2 + 263952c/ - 311040 

°^ ~ 90(c/ + 4)(ci + 2)d(d - 3)(c/ - 7)(d - 9) ' 

(107) 

higher-order coefficients being too cumbersome to write here. The d — > 4 limit of the 
leading coefficient ( |106|) gives the classic result of f2~9fl : |c(l) = — 1| + 0(e). 



Taking the ci — 4 limit of the first 20 coefficients of ( |105| ), we find agreement with p0| , 
where a closed form is given for d = 4. We use these 4-dimensional coefficients to compute 
Pade approximants to the Taylor series. Computation of the two-loop spectral density 
reveals it to be non-negative. Hence II2 is of Stieltjes type and convergence of the Pade 



approximants is guaranteed |23] throughout the region z G [—00, 1] and hence for arbitrar- 
ily large space-like q 2 . Table 8 demonstrates the quality of the convergence. The notation 
is the same as in Section 3.2: truncating the sum after 2n + 1 terms, with the first term 
vanishing, we tabulate E^n), the modulus of the relative error in the 4-dimensional [n/n] 
Pade approximant, for values of q 2 /m 2 up to 12, well outside the disk \q 2 /m 2 \ < 4, to 
which the convergence of the original Taylor series was restricted. 

With just four non- vanishing terms in the Taylor series, the error £^4(2) is less than 
1%, for q 2 /m 2 < 10. The inclusion of each new set of four terms, reduces the error by at 
least two orders magnitude. Tables 6 and 8 clearly demonstrate the utility of applying 
Pade approximants to Taylor series, for obtaining accurate values of two-point functions, 
in the space-like region. Since the calculation of the Taylor coefficients, from vacuum 
diagrams, can be performed in the general mass case IpT] , this is a most welcome result. 



We remark that the technology exists to obtain three-loop QCD results by this method, 
since the calculation of the relevant vacuum diagrams has been automated in the course 
of work on three-loop terms in QED [p6| , [30]. Methods of accelerated convergence may 
also prove useful in the case of multi-loop diagrams with more than two external particles, 
by improving expansions in powers of all the scalar products of the external momenta. 
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Table 8: Accelerated convergence of the Taylor series of II2, with (1 = 4. 



q 2 /m 2 


£4(2) 




£ 4 (4) 




£4(6) 




£4(8) 




12.0 


1.1x10" 


-2 


1.5x10" 


-4 


1.9x10" 


-6 


2.4x10" 


-8 


11.0 


9.4x10" 


-3 


1.1x10" 


-4 


1.1x10" 


-6 


1.2x10" 


-8 


10.0 


7.8x10" 


-3 


7.2x10" 


-5 


6.3x10" 


-7 


5.4x10" 


-9 


9.0 


6.3x10" 


-3 


4.6x10" 


-5 


3.2x10" 


-7 


2.2x10- 


-9 


8.0 


4.9x10" 


-3 


2.7x10" 


-5 


1.4x10" 


-7 


7.6x10- 


-10 


7.0 


3.6x10" 


-3 


1.5x10" 


-5 


5.7x10" 


-8 


2.2x10- 


-10 


6.0 


2.5x10" 


-3 


7.0x10" 


-6 


1.8x10" 


-8 


4.8x10- 


-11 


5.0 


1.6x10" 


-3 


2.8x10" 


-6 


4.5x10" 


-9 


7.3x10- 


-12 


4.0 


8.9x10- 


-4 


8.3x10" 


-7 


7.3x10" 


-10 


6.4x10- 


-13 


3.0 


3.9x10- 


-4 


1.6x10" 


-7 


6.0x10" 


-11 


2.1x10" 


-14 


2.0 


1.1x10" 


-4 


1.2x10" 


-8 


1.3x10" 


-12 


6.2x10" 


-15 



4.4 Accelerated convergence of the asymptotic expansion 

For large q 2 /4m 2 = — we invert the series in ([72|J73| , [74]) and obtain the coefficients of 



n 2 W = E^> ( 108 ) 

n=0 ^ 

for arbitrary d. Each coefficient involves the structures j3 and 7, given in (|9|,|T0|) . In general, 
the terms {1,7,/?, /3 2 } occur, with rational coefficients that are singular as d —>■ 4. The 
d — *■ 4 limit then gives a finite coefficient that is quadratic in ln(g 2 /m 2 ). On the cut, 
with z = —q 2 /4m 2 = t + iO and t G [1, 00], one may extract the asymptotic expansions 
of the real and imaginary parts of n 2 (t + iO), using ln(g 2 /m 2 ) = ln(4t) — m. Here we list 
only the first three asymptotic coefficients, for arbitrary d: 

_8(d 3 - I2d 2 + Aid - M) 32(d - 2)(d 2 - Ad + 8) 

^ ' " ~ d{d - 3)(d - 4)(d - 5) ~ {d-l){d- 4) 3 (3d - 8)(3d - 10) 7 
I 8(d-2)(d 2 -7d+16) 

+ (d-l)(d-3) 2 (d-4)3 /? ' (109) 
16(d-2)(d 3 -8d 2 + 21ci-26) 8(d - 1) 

P ^ ' ~ {d -l)(d- 3){d - 4) 2 {d - 6) (3d - 10) 7 + {d - 3) 2 {d - 4) ^ 
8(d-2)(d 2 -6d+ll)^ 2 ^ 



(d- l)(d-3) 2 (d-4) 
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32(d-l) _ 4(d-2)(d-5)(d 4 -14d 3 + 65d 2 -128d+128) 

P( ^ ~ ~d(d-2)(d-3)(d-4) 2 (d - l)(d - 3)(d - 4) 3 (d - 6)(d - 8) 7 

2(d 3 -3d 2 -26d + 40) n 4(d 3 - 10d 2 + 32d - 36) n2 , , 

^/3H — -3 ill!) 

d(d-3)(d-4) 2 P ^ (d-l)(d-4) 3 P K ' 

Table 9 shows accelerated convergence of the large-g 2 expansion of n 2 , for q 2 /m 2 > 6. 
The error £"4(4) is at least three orders of magnitude smaller than that which would have 
been obtained without accelerating the convergence. As in Section 3.2, the application of 
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the e-algorithm, to both series, allows one to cover the entire space-like region, to good 
accuracy, using only a few terms from each series. Moreover, we find the e-algorithm to 
be similarly effective in accelerating the convergence of the asymptotic expansions of the 
real and imaginary parts on the cut. 

Table 9: Accelerated convergence of the asymptotic expansion of H2, with d = 4. 



q 2 /m 2 


£4(1) 




£ 4 (2) 




£4(3) 




£4(4) 




40.0 


1.1x10" 


-3 


9.7x10- 


-7 


2.7x10- 


-10 


2.3x10- 


-13 


36.0 


1.4x10" 


-3 


1.6x10- 


-6 


5.1x10- 


-10 


5.4x10- 


-13 


32.0 


1.8x10- 


-3 


2.7x10- 


-6 


1.1x10- 


-9 


1.4x10- 


-12 


28.0 


2.5x10- 


-3 


5.0x10- 


-6 


2.4x10- 


-9 


4.1x10- 


-12 


24.0 


3.5x10- 


-3 


9.9x10- 


-6 


5.9x10- 


-9 


1.4x10- 


-11 


20.0 


5.3x10- 


-3 


2.2x10- 


-5 


1.7x10" 


-8 


6.0x10- 


-11 


16.0 


8.6x10- 


-3 


5.8x10- 


-5 


5.4x10- 


-8 


3.3x10- 


-10 


12.0 


1.6x10- 


-2 


1.9x10- 


-4 


1.7x10" 


-7 


2.7x10- 


-9 


10.0 


2.4x10- 


-2 


3.6x10- 


-4 


9.7x10- 


-8 


6.3x10- 


-9 


8.0 


3.9x10- 


-2 


6.5x10- 


-4 


2.8x10- 


-6 


4.9x10- 


-8 


6.0 


7.2x10- 


-2 


1.1x10- 


-4 


1.7x10" 


-4 


2.0x10- 


-7 



5 Summary and conclusions 

There is a pressing need for techniques to calculate multi-loop radiative corrections in 
processes with a wide variety of masses and momenta, in order to probe the standard 
model. Here we have made a second step in harnessing the economy of <i-dimensional 
algebraic methods to that end. The first step was made in || 0, which dealt with 
on-shell massive two-loop two-point functions. Here we have tackled all the off-shell two- 
loop master integrals needed for QED and QCD, using integration by parts to reduce all 
but one of them to hypergeometric functions, and dealing with the troublesome fermion- 
propagator integral of Fig. 1(c) equally effectively, by using differential equations, likewise 
obtained from integration by parts. 

As might have been expected from previous successes 0, [F| §j of dimensional regular- 
ization of massive diagrams, our new d- dimensional results are both easier to obtain and 
simpler in structure than their 4-dimensional specializations. Just as the complexity of 



polylogarithms |22] results from a singular case of hypergeometric series, so 4-dimensional 
calculations are harder than <i-dimensional ones. The moral to be drawn from this is 
clear: take the limit d — > 4 as late as possible. The cost of carrying an extra variable in a 
computer-algebra program is a small one to pay, for the avoidance of analytic complexity, 
and can, in any case, be made even smaller, by retaining only an appropriate number of 
powers of (d — 4). 

We have applied algebraic methods to the two-loop photon self-energy. As promised, 
a thorough knowledge of the analytic properties of the master integral turns out to be 
sufficient to solve problems with arbitrary spinorial or tensorial complexities, without need 
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of further analysis. All problems of the same mass structure (such as the problem ]nj 
of the two-loop gluon-condensate contribution to vacuum polarization, which still lacks a 
full solution) amount to no more than the determination of 5 rational functions of d and 
q 2 /m 2 , multiplying structures that have been thoroughly analyzed in Section 4. 

We have shown examples in which the whole of the space-like region may be accu- 
rately covered by applying Pade approximants to the small-g 2 Taylor series, and their 



generalization, via the e-algorithm | 23fl , to the large-g 2 asymptotic expansion. Tables 8 
and 9 exhibit such accelerated convergence, in the case of two-loop vacuum polarization. 

In the time-like region, our methods are similarly effective. In the case of vacuum 
polarization, accelerated convergence of the Taylor series, for —q 2 G [0,4m 2 ], is achieved 
by applying the e-algorithm to truncations of (|105|) . On the cut, with — q 2 = s + iO 
and s G [4m 2 , oo], accelerated convergence of the asymptotic expansions of the real and 
imaginary parts is achieved by applying the e-algorithm to truncations of ( |108|) . Thus 
we can approximate the two-loop vacuum polarization function (|7S|) to good accuracy, 
throughout the entire complex plane, using only a few of its algebraically computed Taylor 
and asymptotic coefficients. Even near the branchpoint, with — q 2 4m 2 , we find the 
e-algorithm to be effective in accelerating the convergence of the original series. 

These efficient methods of algebraic computation and numerical approximation are 
now available for all the self-energy diagrams of QED and QCD, in all kinematical regions 
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thanks to our explicit large and small q expansions, in (|T7|j20|j24]) and (|T^,1S,22), and our 
analysis of the remaining fermion-propagator integral of Section 2.5, with contributions 
that were conveniently separated in Section 3.1. Given this separation, one may accurately 
approximate the real and imaginary parts in the demanding time-like region — q 2 = s + iO, 
with s G [m 2 ,9m 2 ], by applying the e-algorithm separately to the asymptotic expansion 
of the contribution with the lower threshold and to the Taylor series of the contribution 
with the higher threshold, paying attention to numerical cancellations in the latter's 
coefficients, as in (|61~] , |63D . Elsewhere, one need not make this delicate separation. Tables 6 
and 7 of Section 3.2 show how good is the coverage of the entire region q 2 G [— m 2 , oo], 
thanks to the enlarged domain of convergence of each series. The application of such 



methods to the more general mass cases of [|TT| [12| should help to extend the scope, speed 



and accuracy of standard-model calculations. 

It now requires only careful algebraic programming to reduce the bosonic and fermionic 
two-loop propagators of QED and QCD, in an arbitrary covariant gauge, to the minimal 
set of master integrals of Fig 1, together with a few of their derivatives. The application 
of methods of accelerated convergence to the expansions of these propagators in powers of 
q 2 /m 2 and m 2 /q 2 will then yield economical and accurate parametrizations, throughout 
the entire complex plane. Moreover, terms of order (d — 4) L ~ 2 may also be obtained, for 
use in the dimensional regularization of higher-order calculations, with L loops. 

It seems most likely that some degree of the analytical economy and numerical accuracy 
achieved in this paper is also achievable in wider areas, involving more loops, or more 
masses, or more external particles. We plan to study three-loop diagrams. 
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